Anisotropic magnetic deflagration in single crystals of Gd5Ge4 



S. Velez,i'2^B J- M. Hernandez/- ^ A. Garci'a-Santiago/-^ J. Tejada/-^ V. K. Pecharsky,^.^ 
K. A. Gschneidner, Ji.,^'^ D. L. Schlagel,^ T. A. Lograsso,^ and P. V. Santos^ 

^ Grup de Magnetisme, Departament de Fisica Fonamental, Facultat de Fisica, 
Umversitat de Barcelona, c. Marti i Franques 1, planta 4, edifici nou, 08028 Barcelona, Spam 
^Institut de Nanociencia i Nanotecnologia IN^ UB, 
Universttat de Barcelona, c. Marti t Franques 1, 08028 Barcelona, Spain 
'The Ames Laboratory, U.S. Department of Energy, 
Iowa State University, Ames, Iowa 50011-3020, USA 

^Department of Materials Science and Engineering, 
Iowa State University, Ames, Iowa 50011-2300, USA 
'' Paul-Drude-Institut fiir Festkorperelektronik, Hausvogteiplatz 5-7, 10117 Berlin, Germany 

(Dated: September 7, 2011) 

Experimental evidence of the anisotropy of the magnetic deflagration associated with the low- 
temperature first order antiferromagnetic (AFM) — >■ ferromagnetic (FM) phase-transition in single 
crystals of Gd5Ge4 is reported. The deflagrations have been induced by controlled pulses of surface 
acoustic waves (SAW) allowing us to explore both the magnetic fleld and temperature dependencies 
on the characteristic times of the phenomenon. The study was done using samples with different 
geometries and configurations between the SAW pulses and the direction of the applied magnetic field 
with respect to the three main crystallographic directions of the samples. The effect of temperature 
is nearly negligible, whereas observed strong magnetic field dependence correlates with the magnetic 
anisotropy of the sample. Finally, the role of the SAW pulses in both the ignition and formation of 
the deflagration front was also studied, and we show that the thermal diffusivity of Gd5Ge4 must 
be anisotropic, following hia > > i^c 

PACS numbers: 75.60.Jk, 75.30.Kz, 82.33.Vx 



I. INTRODUCTION 



The Gd5Ge4 intermetallic compound, together with 
other Si doped Gd5(Sia;Gei_2:)4 alloys, has attracted con- 
siderable attention over the last few years, principally due 
to their unusual giant magnetocaloric properties^^"^ . This 
effect is associated with a first-order AFMoFM phase 
transition that occurs simultaneously with a structural 
transformation^'®. The rich phenomenology of these 
transitions in Gd5Ge4 has been broadly studied as a func- 
tion of temperature and magnetic field using both poly- 
crystalline and single crystalline samples2r2^. At temper- 
atures exceeding T 20 K, the isothermal field-driven 
AFM<-)>FM transitions can be continuously reproduced 
but, when the sample is cooled below T ~ 10 K at zero 
magnetic field, the field-driven AFM^>FM transition be- 
comes irreversible. Several reported magnetic behaviors, 
such as the glassy propertiesi^ii^, have suggested that 
the magnetocrystallographic ground state of this system 
at zero magnetic field could be the FM 0(1) [0(1) is a 
common notation of the low- volume orthorhombic poly- 
morph of Gd5Ge4, see Ref. il6i for details], whereas the 
initial AFM 0(11) state [the 0(11) is a common notation 
of the high- volume orthorhombic polymorph of Gd5Ge4], 
which is stable in a zero field cooled (ZFC) process, is 
due to the kinetic arrest of the 0(11) crystal structure. 
However, recent first principles modelingi^r— of the free 
energy of the different 0(1) and 0(H) magnetic phases, 
have pointed out that the ground state of this crystal 
must be the AFM 0(H). 



At low temperatures, the first order field-driven 
AFM— )-FM phase transition can proceed in two different 
fashions. Usually, this transition is rather gradual and 
takes place over a wide range of magnetic fields, but when 
the sample is large enough or the magnetic field is swept 
at high rates, this phase transformation is abrupt, which 
can be identified in the M — H cycle as a magnetic jump^. 
Historically, such magnetic discontinuities have been 
called magnetic avalanches and they have been also ob- 
served in other materiala^Sr— which also exhibit a giant 
magnetocaloric effec t^'^'^^ related to a transition from a 
kinetically-arrested state to the magnetic equilibriuniii^. 
The dynamics of the magnetization of the sample dur- 
ing such transitions have been reported first in molecular 
magnets^"—, and later in manganites^Sr-il and polycrys- 
talline samples of GdsGe^^. For all of these materials, 
it was found that a phase-transition front is formed and 
propagates through the sample at a constant speed on 
the order of a few m/s. The strong similarities between 
this magnetic phenomenon and a chemical combustion^ 
have lead the scientific community to call it magnetic 
deflagration. 

In magnetic deflagration, the role of fuel is played by 
the energy difference between the metastable and the sta- 
ble states of the system, namely AE. This energy differ- 
ence is related to the intrinsic energy of the metastable 
ordered magnetic phase plus the Zeeman energy, that 
comes from the interaction of an external magnetic field 
H with the spins in the system. On the other hand, the 
rate of heat transferred from the region of burning spins 
to their flammable neighbors is controlled by i) the energy 
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barrier to be overcome by the metastable spins, U\ ii) the 
so-cahed characteristic time attempt tq; iii) the thermal 
difFusivity k; and iv) the fraction of flammable spins rim- 
The main difference between chemical combustion and 
magnetic deflagration is that for the latter the source of 
energy is the reordering of the spins of the system instead 
of an irreversible chemical reaction. Therefore, magnetic 
deflagration becomes of special interest due to the non- 
destructive character of the process that allows complex 
phenomena of deflagration to be explored under different 
conditions. 

The theory of magnetic deflagration^^ determines the 
instability condition that leads a typical broad transi- 
tion to the occurrence of a deflagration process. When 
the rate of the thermal jump over the energy barrier for 
a single metastable spin, F, exceeds some critical value, 
Tc, the nucleation of the deflagration front and the sub- 
sequent thermal runaway should take place. This critical 
rate can be written as 



8k{T)kBT^ 



U{H)AE{H)nrnP 



(1) 



where I is some characteristic length, k{T) = k{T)C{T) 
is the thermal conductivity, and C(T) is the specific heat. 
The front of propagation is identified as the flame of the 
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process, whose characteristic size is (5 [K{Tf)/r{Tf)] ' , 
where Tj is the corresponding temperature of the flame, 
which is given by 
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where Od is the Debye temperature. Finally, in the ap- 
proximation of a planar burning-front, the speed of the 
flame is 
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Note that Eq. (1) can be used to establish the thresh- 
old condition for the occurrence of spontaneous magnetic 
deflagration in a field-sweep experimen t'^°i'^'^''^^i^^''^^i'^^ . 
However, in such experiment, it becomes a difficult task 
to explore the field dependence of the speed of the 
flame, v{H). To solve that, experimentalists have de- 
veloped techniques to study the magnetic deflagration 
using controlled experimental conditions. These experi- 
ments consist, basically, of sending a heat pulse that acts 
as a spark of flame that ignites the deflagration. At- 
tached resistors^2iMi2^i^^ electrical contacts made on the 
sample^ or surface acoustic waves (SAW)^i^ are exam- 
ples of sources that can be used for this purpose. 

However, the test of magnetic deflagration has been 
limited to only a single law of propagation. In the case 
of molecular magnets, the speed of the deflagration front 
is determined by the value of the magnetic field applied 
along the easy magnetization axis, whereas the trans- 
verse field affects the threshold conditions [see for exam- 
ple ref. [13 via their unusual quantum propertie a^^'^^ . In 



the case of manganites, as well as polycrystalline samples 
of Gd5Ge4, there is no influence, excluding geometrical 
effects, of the direction of the applied magnetic field on 
the properties of the defiagration process, whose observed 
characteristics are the result of averaging the properties 
along the principal crystallographic axes of the sample 
due to their random distribution. Therefore, the goal of 
this work is to investigate whether the magnetic defla- 
gration in single crystals of Gd5Ge4 is anisotropic, and 
what is the role of each crystallographic axis in both the 
formation and propagation of the deflagration front. 



II. EXPERIMENTAL SET-UP 

A large (approximately 4 mm diameter, 40 mm long) 
single crystal of Gd5Ge4 was grown using the tri-arc 
technique^. The Gd metal used to prepare the sto- 
ichiometric polycrystalline charge weighing 20 g total 
was prepared by the Materials Preparation Center of the 
Ames Laboratory^, and it was at least 99.99 wt.% pure 
with respect to all other elements in the periodic table. 
The Ge was purchased from Meldform Metals, and it 
was 99.999 wt.% pure. The as-grown single crystal was 
oriented using backscatter Laue technique. Two different 
single crystals of Gd5Ge4 used in this work were cut from 
a larger single crystal using spark erosion. Their dimen- 
sions were as follows: sample 1: Si « 1.17 x 2.45 x 1.04 
mm'^ and sample 2: 5*2 ~ 2.40 x 1.29 x 1.07 mm^ along 
the crystallographic directions a, 5, and c, respectively. 

Figure 1 shows the schematic of the experimental set- 
up used in our measurements. The Gd5Ge4 samples were 
mounted using non-magnetic commercial silicon grease 
on a piezoelectric device [Fig. 1(a)] specially designed to 
send SAW pulses to the sample. The excitation spectrum 
of SAW modes with this system is basically determined 
by the resonances of the interdigital transducer (IDT), 
whose values are multiple harmonics of a fundamental 
frequency /o « 111 MHz (see, for example, ref. 21,, 
and m for more details). The device is placed inside 
a commercial SQUID magnetometer [Fig. 1(b)] capa- 
ble of measurements at temperatures down to 1.8 K in 
magnetic flelds up to 5T. The microwaves, for the SAW 
generation, were transported from an external commer- 
cial Agilent signal generator to the IDT placed inside 
the cryostat by means of a coaxial wire that introduces 
attenuation smaller than 10 dB. 

The SAW pulse induced in the IDT propagates along 
the lenght of the piezoelectric substrate made of LiNbOs. 
The crystallinity of the substrate, together with the ge- 
ometry of the device, provides an amplitude profile of the 
SAW oscillations in the XZ plane, being therefore the 
direction of oscillation of the SAW waves out of plane, 
that is parallel to the y direction [see Fig. 1(a) for the 
definition of the axes] . Since z is the largest side of the 
LiNbOa crystal and x is the perpendicular direction, the 
amplitude profile of the SAW can be considered prac- 
tically independent of z, but it depends on x having a 
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maximum at a; = (that is, t he p osition of the center of 
the IDT; see for example ref. [50, and references therein, 
for more details). When desirable experimental condi- 
tions (a specific combination of T and H) are reached 
and stable, a controlled SAW pulse is delivered to ignite 
the magnetic deflagration process in the Gd5Ge4 sample. 



(a) 1 
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FIG. 1: Experimental set-up. (a) Schematic of the piezoelet- 
ric device used. 1, coaxial cable; 2, conducting stripes; 3, IDT; 
4, LiNbOs; 5, single crystalline sample of Gd5Ge4. (b) View 
of the spatial distribution of the piezoelectric set-up inside 
the SQUID magnetometer. 1, coaxial cable; 2, piezoelectric 
set-up; 3, sample holder; 4, superconducting coils for the field 
generation; 5, magnetometer's pick-up coils. 

The magnetic field was always applied along the z di- 
rection of the sample holder [as shown in Fig. 1(b)], and 
the piezoelectric set-up was placed between the pick-up 
coils of the SQUID magnetometer. Two different tech- 
niques have been used to obtain magnetic measurements 
with this system. The first, and the typical mode, con- 
sists of taking dc magnetic measurements by moving the 
sample through the four pick-up coils of the magnetome- 
ter. It allows one to obtain the absolute magnetic mo- 
ment of a given sample with a very high precision. This 
technique has been used i) to characterize the magnetic 
properties of the sample and ii) to verify the magnetic 
state of the sample before and after each induced defla- 
gration. The limitation of this mode is the time it takes to 
perform each measurement, which is approximately ~ 30 
seconds. The second method consists of directly mea- 
suring the voltage from the SQUID-voltmeter without 
requiring any motion of the sample. Placing the center 
of the sample at the position of the inner coils [as shown 
in Fig. 1(b)], where the sensitivity of the system has its 
maximum, the voltage drop recorded is directly related to 
the magnetic state of the sample [i.e., AI/ oc AAf]. This 



technique allows us to monitor fast magnetic changes 
with a time-resolution better than 0.01 ms, so that we 
can measure the time evolution of the magnetization of 
the sample during a magnetic defiagration process. 

Each single crystal has been studied under different 
sample/set- up configurations to elucidate the role of each 
crystallographic axis on the properties of the defiagration 
phenomenon. Due to geometrical restrictions, four differ- 
ent configurations were available to be explored for each 
sample. When the magnetic field is applied along the 
longest side of the crystal, the SAW pulse can be applied 
to either of the two shorter sides. On the other hand, if 
the applied field is along one of the shorter directions, the 
SAW pulse can be applied parallel to the other short side. 
We will refer to each sample configuration using the fol- 
lowing notation: Si{x,y, z), where i denotes the sample 
number, and x, y and z, correspond to the orientation of 
each crystallographic axis of the sample with respect to 
the coordinate system of the sample holder shown in Fig. 
1(a). For example, 82(0, b,c) refers to sample 2 with its 
a-axis parallel to the cc-direction of the sample holder, b- 
axis parallel to the y-direction, and c-axis parallel to the 
z-direction, which is the direction of the magnetic field 
vector. 



III. RESULTS 

DC magnetic measurements of the field-driven 
AFM— i>FM transition at low temperatures were carried 
out for both Gd5Ge4 samples and for each configuration 
after the sample was ZFC from T = 50 K. Figure 2 shows 
the AFM— J'FM transformation and the subsequent re- 
moval of the magnetic field in the FM state obtained 
at r = 2 K for each independent crystallographic axis 
of the system using sample 1. The same results were 
obtained for sample 2 but these are not shown here for 
brevity. As the magnetic field is increased from zero, 
the linear slope of the M{H) curves [note that the ob- 
served step around H ^ 8 kOe when the magnetic field 
is applied along the c-axis is associated with the spin-fiop 
transition^] suggest that the initial state of the sample is 
purely AFM, and it remains unchanged, until direction- 
specific critical magnetic field. He, is reached, whose val- 
ues at T = 2 K are 28 kOe, 23 kOe and 26 kOe for the 
a, b and c crystallographic axis, respectively. Above it, 
the AFM— J^FM transformation is quite gradual, and it 
takes place over a field range AH ^ 4 kOe. The inset 
of Fig. 2 shows, specifically, the fraction UAFAiiH) of 
metastable AFM spins around the transition. At higher 
temperatures, the critical field Hc{T) decreases with the 
rate of dHrJdT « -1.5 kOe/K in the range 2 - 8 K. All 
of these results are in agreement with the previous data 
reported for single crystalline Gd5Ge4 samples^i. 

After each crystallographic axis has been magnetically 
characterized, we proceeded to perform the deflagration 
experiments. Two control variables were i) the initial 
temperature Tq and ii) the applied magnetic field for the 
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to induce transformation of only a fraction of spins into 
the stable FM state without deflagration. Moreover, the 
fraction of transformed spins decreases with decreasing 
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FIG. 2: (color online) Magnetic anisotropy of the field-driven 
AFM— >FM transition and the subsequent removal of the ap- 
plied magnetic field measured at To = 2 K after sample 1 was 
ZFC from T = 50 K. The inset shows the fraction uafm of the 
metastable AFM spins around the magneto-crystallographic 
transition. 
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FIG. 3: (color online) Normalized time evolution of the 
change in the magnetization of the sample, AM (t) / AMt , in 
the case of Si (c, a, b) for different ignition fields Hig at T = 2 
K at which a magnetic deflagration was observed. From right 
(red) to left (blue), the signals shown correspond to the data 
obtained for the values explored from Hig = 16.5 kOe to 
Hia = 22.0 kOe, respectively, in steps of 0.5 kOe. For each 



ignition Hig. For each deflagration measurement, the 
sample was ZFC from T = 50 K to the desired tem- 
perature Tq. After that, the applied magnetic field was 

increased slowly to a selected ignition field H,g < H^ and measurement, the sample was first ZFC from T = 50 K, and 
then a SAW pulse of 100 ms width and 16 dB was deliv- then a H,g was apphed. After that, a SAW pulse of 100 ms 
ered. As a consequence of the pulse, the sample is driven width was delivered at t = 0. 
out of the initial equilibrium and for a certain range of 
experimental Tq and Hig conditions, a magnetic deflagra- 
tion was induced. The time evolution of the change in 
the magnetization of the sample, AM{t), was recorded 
from the SQUID-voltmeter, where t — corresponds to 
the delivery of SAW pulse. 

Figure 3 shows the resulting AM (t) recorded for 
Si{c,a,b) configuration at To = 2 K in the range of ig- 
nition fields at which the occurrence of magnetic defla- 
gration was identified, and the initial value of uaf m was 
kept close to one. The data were normalized to the total 
magnetic drop observed for each case, AMt, which corre- 
sponds to the variation of M when the spins of the sample 
change from the AFM state to the full FM state. This 
was confirmed from dc-magnetic measurements taken be- 
fore and after each SAW pulse was delivered. For lower 
ignition fields [Hig < 16.5 kOe], the same kind of mea- 
surements performed revealed that no more than ^ 10% 
of the spins of the sample become FM (not shown for 
simplicity). This abrupt change of the number of spins 
that transforms to the FM state with Hig is a typical 
feature that shows the self-maintenance of the deflagra- 
tion process that utilizes the energy of the metastable 
spins for the occurrence of the full phase transformation 
when the threshold for the ignition of the deflagration is 
reached. In other words, in our experiment this threshold 
is reached for ignition fields from 16.5 kOe to 22.0 kOe, 
but at any smaller field, the effect of the SAW pulse is 



From the data shown in Fig. 3, two characteristic times 
related to the deflagration phenomenon can be identified 
for each curve: tig, that is the time that is required to 
reach the deflagration threshold after the SAW pulse has 
been switched on, and corresponds to the point at which 
the slope in the AM(t) curve exhibits an abrupt change; 
and td, deflned as the subsequent time interval at which 
this fast change takes place due to the magnetic defla- 
gration. Note that both characteristic times are strongly 
influenced by the ignition fleld, but no remarkable dif- 
ferences were observed when the experimental procedure 
was repeated at fixed Hig for a wide range of different Tq 
values (not shown). It is an expected behavior that fol- 
lows from the theoretical expressions given above. From 
Eq. 3, v{H), which is related to td{H), only depends 
on H and is independent on Tq, on condition that rim 
is constant. This is our case because we restrict our ex- 
periments to Um = nAFM — 1- On the other hand, the 
value of tig can be theoretically found solving the inequal- 
ity T{H,T) > Tc{H,T). Essentially, it is accomphshed 
when T is increased enough due to the delivery of a SAW 
pulse (in section IV. B we will show why the effect of the 
SAW pulse is to heat the sample), providing a certain tig 
value for each experimental condition. However, testing 
this equation for different Tq and Hig it can be shown 
that, whereas a small change in Hig provides a strong 
change in tig , the effect of Tq is practically negligible for 
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FIG. 4: (color online) The times of deflagration, td, and the ignition times, tig, obtained at different Hig for the two samples 
studied under different sample/set-up configurations split in two groups, (a) and (b) corresponds to td and tig, respectively, 
when the applied field is along the longest side of the sample, (c) and (d) corresponds to the same data taken but for the 
perpendicular configurations. 



a wide range of experimental values. Additionally, the 
higher Tq , the smaller the capability of the system to re- 
tain the spins in the metaestable AFM state, reducing 
the experimental range over vifhich the phenomenon of 
magnetic deflagration can be explored. Considering this, 
we have focused on the field dependencies of the proper- 
ties of the deflagration phenomenon among the different 
crystallographic axes and sample configurations. 

Figure 4 shows the characteristic times, td{Hig) and 
tig{Hig), obtained for the two samples and for the differ- 
ent configurations explored split in two groups accord- 
ing to different geometrical arrangements of the crystals 
in the sample holder: magnetic field applied along the 
longest side of the crystal [Fig. 4(a) and 4(b)] and mag- 
netic field applied perpendicularly [Fig. 4(c) and 4(d)]. 
To avoid both the ignition and the extinction of the de- 
flagration process that do not correspond well with the 
deflagration of a flame through the sample, we have used 
the middle part of the total magnetic change to identify 
td- In particular, we have taken the time elapsed between 



25% and 75% of AMt multiplied by 2 to rescale the 
data. Independently of the specific characteristic times 
observed for each configuration, several common features 
are worth noting. When Hig is close to He, related to the 
crystallographic axis of the sample that is parallel to the 
applied magnetic field, both the ignition and deflagra- 
tion times are rather fast taking place only in a few tens 
of ms and few ms, respectively. However, as Hig moves 
down and away from He, both times increase progres- 
sively showing a non-linear dependence. Moreover, when 
Hig is reduced far from He, the energy supplied by the 
SAW pulse is no longer sufficient to reach the ignition 
threshold, and therefore, the deflagration does not take 
place. Moreover, depending on the configuration of the 
sample, a few deflagrations take place even after the SAW 
pulse has been switched off [for example, see Si{b,a,c) 
in Fig. 4(d) where five deflagrations exhibit tig > 100 
ms], indicating again the self-maintenance character of 
the deflagration phenomenon. 
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IV. DISCUSSION 

There are two important features to take into account 
from the time values shown in Fig. 4. On one hand, 
for a certain ignition magnetic field Hig, the observed 
values of both characteristic times of the magnetic defla- 
gration mainly depend on the crystallographic direction 
of the sample along which the magnetic field is applied. 
For example, this follows from strong differences observed 
between the values obtained for a and b axes, where the 
field intervals at which the deflagrations take place do not 
overlap. However, notice that if the relative field value 
Hig — He, where He is the corresponding critical field 
along which the magnetic field is applied, is taken into 
account, the characteristic times observed become similar 
indicating that there is a correlation between the mag- 
netic anisotropy of the sample and the observed deflagra- 
tion times. On the other hand, the differences between 
the sets of data Si{a,c,b) and Si{c,a,b) [or £'2(6,0,0) 
and 82(0, b, a)], where the geometrical dispositions of the 
sample are equivalent (which would imply same lengths 
of propagation) and the applied magnetic field is along 
the same crystallographic axis, it suggests that additional 
effects on the properties of the magnetic deflagration may 
come from the crystallographic orientation of the sample 
in the XY plane of the sample holder. Therefore, the 
aim of this section is to discuss these experimental obser- 
vations within the framework of the theory of magnetic 
deflagration to show i) the connection between the mag- 
netic anisotropy of the sample and the observed field de- 
pendencies and ii) the role of each crystallographic axis 
of the sample and the SAW pulse on the ignition and 
propagation of the flame. 



A. Comparison of the data with the Theory of 
Magnetic Deflagration 

For simplicity, and because td{H) is a more reliable 
fingerprint than tig of the deflagration phenomenon, we 
start our discussion focusing on these values. They must 
correspond to td{H) ~ lp/v{H), where we have deflned 
Ip as the length along which the deflagration front prop- 
agates inside the sample and v{H) is its speed. For a 
planar front propagating along one of the principal crys- 
tallographic axis of the sample, the time of propagation 
of the flame must follow: 




UjH) 
Tf{H) 



exp 



UiH) 



kBTf{H) 



(4) 



Omitting nonessential factors that cannot contribute 
to any observable difference of the values of td{H) be- 
tween the samples and configurations, all the experimen- 
tal dependencies to be taken into account can be split 
in two different functions according to their origin: the 



term f{H) is related to the magnetic field dependencies, 
i.e., the energy barrier U{H) and the temperature of the 
fiame Tf{H) [see Eq. (2)], which essentially depends on 
AE{Hy/^ because in all of the data reported tiafm — 1- 
Notice that these field-dependent values correspond to 
the crystallographic axis along which the magnetic field 
is applied. On the other hand, geometrical contributions 
associated with the heat transport have been grouped in 
g{lp, k)'. the propagation length Ip and the thermal diffu- 
sivity K, whose values correspond to the crystallographic 
direction along which the flame propagates. In the case 
that it does not propagate only along one direction, the 
g{lp, k) function should be modified with an appropri- 
ate combination of the values of k and Ip for the axes 
involved in the process (we will further discuss this ques- 
tion in section IV. B). The characteristic time attempt. 
To, has been considered as a constant because it is typi- 
cally a global parameter that characterizes the time-flip 
of all the spins of a solid and it should be independent of 
the direction at which the magnetic field is applied. 

For any set of data, the main contribution to td{H) 
comes from the f{H) function due to the strong field- 
dependent term exp[C/(iJ)/fcs77 (iJ)]. Notice that this 
exponential dependence can explain the non-linearity ob- 
served in Fig. 4. Since Tf{H) oc A£'(iJ)i/'', the values 
of td{H) should be basically determined by the field de- 
pendence of the energy barrier, U{H), whereas Tf{H) 
may be considered as a constant for each set of data due 
to the limited range of experimental fields explored. Ad- 
ditionally, due to the weak magnetic anisotropy of the 
material, the values of Tf should be similar for all crys- 
tallographic directions in our experiments. In this con- 
text, both the scaling of td through the relative value of 
Hig with respect to the corresponding He, and the mag- 
netic anisotropy of the sample, can be directly explained 
through the anisotropy of the field dependence of the en- 
ergy barrier, U{H). To better understand it, for a given 
temperature and crystallographic direction, notice that 
He correspond to the field at which the effective energy 
barrier U {H) , is reduced enough to get the characteristic 
time attempt r = tq exp[?7(if)/fcBT] of the order of the 
time of measurement^. According to this, and taking 
into account that the field dependence of U should be 
the same for all crystallographic directions, for same rel- 
ative field values, Hig — He, the strength of the energy 
barriers would be equal for the three axes, and therefore, 
similar values of td are expected despite the differences 
that could arise from the geometrical term. 

To be more precise in our further discussion on the geo- 
metrical contribution, g{lp, n), to the defiagration values, 
we proceed to compare the data obtained for different 
configurations when the applied magnetic field is along 
the same crystallographic direction. In this case, the sets 
of data are controlled by the same f{H) function. Ad- 
ditionally, for a similar geometrical dispositions of the 
sample in the sample holder, Ip is expected to be similar. 
Therefore, the observed differences [for example, see Fig. 
4(a)] must be attributed to an anisotropy of k. However, 



7 



to better understand this fact, we have to determine, 
first, what is the role of both the SAW pulse and sam- 
ple configuration on the nucleation of the deflagration 
front and, second, through which length(s)/direction(s) 
the deflagration front propagates. It could also give a rea- 
son for the different observed tig{Hig) values. However, 
before that, it is interesting to know whether the predic- 
tions of either the speed of the deflagration front or the 
values of the ignition time that arise from the theory of 
magnetic deflagration^ are in the range of the observed 
experimental values. 

From previously reported datap^ of the magnetic de- 
flagration in a polycrystalline sample of Gd5Ge4, Tf was 
estimated to be around 30 K. Taking into account the 
discussion given above about the Tf{H) values of the dif- 
ferent crystallographic axis, and that the magnetic prop- 
erties of a polycrystal are averaged, we can assume that 
T/ ^ 30 K is a good estimate of the temperature of the 
flame in all observed deflagrations. Thus, the thermal dif- 
fusivity of the sample is estimated to be^ i^{Tf) ~ 3-10~^ 
m^/s. Taking /p ~ 1 mm, tq ~ 10~^ s, and that in our ex- 
perimental field range U (H) should be around 200 — 300 
K, one gets from Eq. (4), td ^ 10 ms, which is in a 
good agreement with our observations. The other in- 
teresting parameter to test is the width of the fiame S, 

1/2 

which can be found through 5 - [K{Tf)/r{Tf)Y'. Here, 
T{Tf) = TQ^exp[—U/kBTf], so the upper limit for the 
width of the flame is estimated to be Smax 0.1 mm, 
which means that such a flame can be formed and can 
propagate inside the sample. 

For the ignition time values, they can be estimated 
solving the inequality T > Fc [Eq. (1)]. From Ref. HH, 
fc(T) ^ 8 J/s-K. Using all of the estimated values given 
above, the condition F ^ Fc should be accomplished for 
Tig ^ 12 K. In other words, the deflagration process may 
take place if the temperature in some part of the sam- 
ple is quickly raised above this threshold temperature^. 
Although in our experiments this cannot be done fast 
enough, we may assume that this condition could be ac- 
complished if a certain volume of spins of the sample, of 
the order of the spins contained in the burning front, are 
heated around this temperature. Since the upper limit 
for the width of the deflagration front is 6max 0.1 mm, 
the upper limit for this volume is on the order of 10% of 
the total volume of the sample. Then, we proceed to esti- 
mate the thermal rise experimented by this volume due to 
a SAW pulse of At = 100 ms. Assuming a good transfer 
between the SAW pulse and the sample, the transferred 
energy is estimated to be Esaw = P ■ St ^ 10 mW • 0.1 
s = 1 mJ. On the other hand, the heating of Ng mols can 
be expressed as dT — dEsAW /C{T) ■ Ng. In our case, 
the mass to be heated is ms ^ 1 mg, so that Ns ~ 10~® 
mols. In the region of interest, C{T) = aT^, where a 
is a well-known constant whose value is^ 0.7 J/mol-K'*. 
Therefore, the final temperature of this volume is given 
by Tg = AAEsAW /oim-s + Tg, where Tq is the initial 
temperature. Replacing all parameters in this equation 
with the known values, one gets ~ 9 K, which is in 



agreement with the expected ignition temperature. Tig. 



B. Formation of the deflagration front 

Due to the location of the sample with respect to the 
piezoelectric device, and taking into account both the 
geometry of the sample and the quasi independence of 
the amplitude of the SAW oscillations in the z direction, 
the spatial description of any deflagration phenomenon 
should be described in two dimensions, in which both 
the ignition and the propagation of the deflagration can 
be described in the XY plane. Taking into account the 
profile of the SAW oscillations in the x direction of the 
piezoelectric device, it is reasonable to assume that the 
energy is supplied to the sample mainly at the center of 
the bottom surface, defined as the (0,0) point in the XY 
plane. On the other hand, at the range of our working fre- 
quencies, the phonon thermalization process should oc- 
cur in less than 1 ms, whereas the characteristic times 
involved in our experiments are at least of the order of a 
few ms. All these features suggest that the interaction of 
the SAW pulses can be approximated as a spark of fire 
that essentially heats the sample at the (0,0) point. 

The heat supplied to the sample, during a SAW pulse, 
diffuses in the XY plane resulting in a thermal rise that 
depends on the position and the time elapsed from the 
ignition of the pulse. Essentially, each isotherm follows 
an ellipsoidal shape in the XY plane, whose characteris- 
tic lengths should follow the relation l^/ly ~ {^x/ HyY^"^ , 
where k denotes the length of the i axis of the sample re- 
ferred to the (0,0) point and Ki is the thermal diffusivity 
of this axis. In the case of isotropic diffusion, the char- 
acteristic lengths Ix and ly should be equal. Therefore, 
any deflagration ignited in any orientation of the sample 
in the XY plane has to exhibit the same characteristic 
times. However, if one crystallographic orientation ex- 
hibits higher thermal diffusivity, the supplied heat pen- 
etrates easier in such direction breaking the symmetry 
of the plane and, consequently, the resulting deflagration 
properties. 

To illustrate the effect of the anisotropy of the thermal 
diffusion in our experiments, let us show what should 
happen with the characteristic times of the magnetic de- 
flagration in the case of a square geometry [L^ = Ly, 
which corresponds to the case of the data shown in Fig. 
4(a) and Fig. 4(b)] but with a strong anisotropy of the 
thermal diffusion [for example, suppose ki >> K2]. It 
can be easily found that the difference between the igni- 
tion times, when ni is oriented parallel to the x direction 
in comparison with the resulting times when it is oriented 
along the y direction, approaches tig{K^y) /tig{Ki\^x) 
2 because, whereas in the first case the phase front should 
be generated when Lx/2, in the second one ly — t- Ly 
due to the point at which the SAW pulse is transferred 
to the sample. Therefore, the distance to be covered by 
the defiagration front formed is different in each case, 
with the ratio between the deflagration times approach- 
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ing td{ni\\y)ltd{ni\\x) 1/2. 

Focusing on the data shown in Fig. 4(a), when the y 
direction is along the c crystallographic axis [Si (a, c, h) 
and S2{h, c, a)], td is higher than when either 5 or a axis 
are ahgned in this direction [6*1(0,0,6) and 62(0, 6, a)]. 
From the phenomenological point of view, the shape of 
the td curves are the same when the y direction is par- 
allel to c, and on the other hand, the curves at which 
c is parallel to x are also similar between them. More- 
over, Fig. 4(b) shows that for these sets of data, higher 
tig values are obtained when c is parallel to y direction. 
These findings, together with the expected thermal dif- 
fusivity dependencies previously discussed, are indicative 
of smaller thermal diffusivity along the c axis compared 
to b and a axis. From similar arguments, the values ob- 
tained in other configurations [for example, the data cor- 
responding to the sets Si{a,h,c) and 62(6, a, c)], suggest 
that Ka should be higher than k^. 

Finally, note that, for H applied along the same crys- 
tallographic direction, the differences of td{H) values ob- 
served between configurations cannot be explained by a 
unique scaling factor. It should be attributed to a dif- 
ferent ratio of distances covered by the deflagration front 
depending on the ignition field explored. At magnetic 
fields close to He , the system can be driven out of equi- 
librium easier. In other words, the deflagration front 
boundary should form close to the (0,0) point, whereas 
for smaller fields, its formation is more complicated and 
it should occur deeper inside the sample where the igni- 
tion can take place^i^. Therefore, when H ^ He, the 
phase front propagates all over the XY surface. However, 
when H explored is far away from He, the anisotropy 
of K implies that the defiagration front formed should 
be different depending on the configuration studied as 
has been explained before, and then, different propaga- 
tion lengths and deflagration times are expected between 
them. In conclusion, the different observed fleld depen- 
dencies should be attributed to the geometric function, 
g{lp, k), as a consequence of non-trivial interplay between 
how the front is generated and how it diffuses in the XY 
plane. 

V. CONCLUSIONS 

To summarize, magnetic deflagrations associated with 
the first order AFM— >-FM magneto-crystallographic 
transformations in single crystals of Gd5Ge4 have been 
induced by controlled SAW pulses. The study has been 
done for different experimental conditions and configura- 
tions between the SAW pulses and the applied magnetic 
field with respect to the crystallographic axes of the sam- 
ples. As expected, the dynamics of the process fits well 
within the framework of magnetic defiagration theory, 
but the comparison of the data obtained between differ- 



ent configurations have revealed anisotropic character of 
the process associated with both magnetic and thermal 
properties of each of the three crystallographic axes of 
the sample. The main effect comes from the field depen- 
dence, which is correlated with the magnetic anisotropy 
of Gd5Ge4 through the anisotropic character of the field 
dependence of the energy barrier, U{H). 

The data obtained suggest that the thermal diffusiv- 
ity is anisotropic, following Ka > Hb > i^c- It plays 
an important role in the front formation and the sub- 
sequent propagation inside the sample due to the fact 
that the anisotropy of the thermal diffusion can be in- 
terpreted as hard and/or easy axes for the occurrence of 
the phenomenon. Reported previously anisotropy of elec- 
trical resistivity and magnetoresistance^, sound propa- 
gation and elastic properties'^ in these and related al- 
loys, makes the conclusion about anisotropy of the ther- 
mal conductivity reasonable. However, additional exper- 
iments should be performed to confirm and further ex- 
plore this property. 

The role of the SAW pulses in the ignition of the mag- 
netic defiagration has been also highlighted. The direc- 
tionality of the SAW pulse transferred to the sample and 
the characteristic times evolved in the deflagration pro- 
cess suggest that the SAW pulses act as an unidirectional 
heater leading to the deflagration process ocurring in the 
perpendicular cross section of the sample. However, we 
note that for systems in which tig{H) or td{H) are less, 
or at least, of the order of 1 ms, the phonon-spin interac- 
tions could play an important role in the properties of the 
magnetic deflagration^. Finally, while this work concen- 
trates on the anisotropy of the dynamics of the deflagra- 
tion phenomena, the authors want to remark that simul- 
taneous to the magnetic deflagration process, a structural 
change takes place in the system. Further studies should 
elucidate the very interesting physics of what is happen- 
ing, principally, inside the magneto-structural burning 
front. 
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